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In the framework of the recently proposed saturation picture, we examine in a systematic way 
whether the nuclear modification factor measured for d-Au collisions at RHIC may be simply ex- 
plained. The Cronin peak which is obtained at mid-rapidity around k± ~ 3 GeV may be reproduced 
O |' at the proper height only by boosting the saturation momentum by an additional nuclear component 

as already shown in the literature. In this respect, mid-rapidity RHIC data cannot necessarily be 
seen as a probe of the saturation picture. The large rapidity (77 ~ 3) region allows us to test the 
shape of the unintegrated gluon distribution in the nucleus, investigating various parameterizations 
inspired by large rapidity solutions (of the BFKL and) of the Balitsky-Kovchegov (BK) equation. 
A satisfactory description of Rcp at RHIC is obtained in the BK picture. 

X' 



Abstract 



c5 ' 1 Introduction 



Testing the saturation (Color Glass Condensate) picture pp, for the initial state of deuteron-gold (dA) 
collisions against RHIC data has been a subject of interest for some time. Two salient features have 
been observed 0E1 concerning the behavior of the nuclear modification factor RdA- At mid-rapidity, 
the Cronin peak height depends on the centrality of the collision. At large rapidity, the suppression 
predicted by quantum evolution is observed and is bigger for smaller centralities. A number of papers 
have recently discussed the description of the Cronin enhancement ^ 03 El EJ |5J |U] and the effect of 
small x HO El I1R 1121 1131 HI) . In this work, we examine in a systematic way, how and if saturation and 
quantum evolution provide a reasonable quantitative agreement with data |2]- 

The time is indeed appropriate to assess the predictability of the saturation (CGC) picture. This 
endeavor has a mitigated conclusion: as it turns, unavoidably, the saturation scale introduced in the 
theory, does not have the proper size to explain RHIC data at mid-rapidity. This conclusion is similar 
to the one stated in [S]. On the other hand, quantum evolution as described by the theory gives the 
proper suppression of the nuclear modification factor above the saturation scale. 

* E-mail address: mehtar@th.u-psud.fr 
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In section 2.1, we calculate the hadron production cross section in dA at mid-rapidity using the 
semi-classical approach and show the prediction for the minimum-bias nuclear modification factor RdA 
(and Rcp for central versus peripheral collisions dA collisions) in relation with other previous studies 
and comparing with data. 

We then discuss, in section 2.2, quantum evolution. We first derive the expression for the cross-section at 
leading log accuracy, including both gluon and quark distributions within the deuteron. This expression 
is identical to eq. (22) in (15) . We then present various parameterizations of the unintegrated gluon 
distribution in the nucleus, inspired by large rapidity solutions of the BK equation (2J and show the 
comparison with representative data j2j. In section 3, the conclusion and outlook are given. 



2 Hadron production in dA 

The nuclear modification factor RdA and the Rcp (Central/Peripheral collisions) ratio are defined as 



j dA — > hX 



dN a 

dA AT ,, dNPP^ hx ' > 

p dN dA^hX 
"coll Airf^k C 



j-, coll drid' 2 k i~ /(~.\ 

Rcp - nC W 



oil dijd 2 k \r 

k and r\ are respectively the transverse momentum and the pseudo-rapidity of the observed hadron. N co u 
is the number of collisions in dA, it is roughly twice the number of collisions in pA(proton-Gold). The 
centrality dependence of RdA is related to the dependence of j\f dA ^ hx = da dA ^ hx / d 2 h and N co u(b) 
on the impact parameter of the collision. In this paper, we address the predictions of the Color Glass 
Condensate for these ratios. We always assume that cross-sections depend on the impact parameter 
only through the number of participants which is proportional to the saturation scale 



} 2 sA 



(b) ~ Q 2 sA {0)N part . Au {h)/N p 

art.Au (0), (3) 



where N part .Au is the number of participants in the gold nucleus in d-Au collisions. This is coherent 
with the assumption that Q 2 sA (b) ~ (N part .Au(b) /2)Q 2 sp such that Q 2 sA (b = 0) ~ A 1 / 3 Q 2 sp [T7]. We use 
Table 2 in [2] which gives the number of participants N par t and the number of collisions N co u for several 
centralities. 



2.1 Semi-classical approach 

We first deal with gluon production at mid-rapidity for which different approaches have been proposed. 
The inclusive cross section has been calculated in ^S] in a quasi classical approach of multiple rescattering 
inside the nucleus, see also H3EI1- Further confirmation has been given in The inclusive cross 

section, for a gluon with transverse momentum k and rapidity r/ — 0, is written as 

rf—dA—^gX (~i i r r 

S^PEPE - J £B J «i^V> (z,b-B)V= Wo (z,b )e -«", (4) 

where Nq(z, b) is the forward scattering amplitude of a gluon dipole off the nucleus, it contains all 
higher twists in the semi classical approximation and i%q(z, b — B) is the forward scattering amplitude 
of a gluon dipole off the deuteron at leading twist approximation. B and b are the impact parameters 
of the deuteron and the gluon with respect to the center of the nucleus. At mid-rapidity and RHIC 
energies, it is legitimate to neglect quantum evolution. 
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This approach provides a description of N G at all orders in terms of the saturation scale Q 8 a ~ 
Q 2 cc A 1 ' 3 Aq Cd . The CGC approach yields a Glauber-Mueller form for the dipole forward scattering 
amplitude 

iV G (z,b) = l-exp(-iz 2 Q2 ( b)ln-i^), (5) 
where the saturation scale is given by [T%] 

1 47r 2 rT 

Q 2 (b) In = -^ pT (b)xG(x, 1/z 2 ). (6) 

p is the nuclear density in the nucleus and T(b) = 2\J R 2 — b 2 the nuclear profile function of a spherical 
nucleus of radius R. A is an infrared cut-off of order Aqcd- The cross section can be rewritten in the 
following simpler form 

t^k = q ^il ' duln^d u [ U d u N G (u,b)}MMu), (7) 

where u = |z|. 



2.1.1 A model for minimum bias collisions 

One defines the minimum bias cross-section as the average over the impact parameter of the collision, 
it may be written as 

dN min.bias d(J 1 f o d<T 



dr]d 2 k x drjd 2 kd 2 b' S A J ^ b drjd 2 kd 2 b' ^ 

Sa = kR 2 is the transverse area of the nucleus. All the centrality dependence in the cross-section is 
contained in Ng(z, b) as given by (JSJ) and I©- One can perform the b integral for Nq yielding 



(N G (z, Q,)) = 1 + — [(1 + -z 2 Q 2 sC ln(l/z 2 A 2 )) exp{--z 2 g 2 c ln(l/z 2 A 2 )} - 1], (9) 

z QsC^ n (l/ z A z ) a o 



where Q 2 C = Q 2 (b = 0). The corresponding result for RdA is shown in Fig. (a) taking Q 2 C = 2 
GeV 2 and A = 0.2 GcV. The proton-proton cross-section is calculated by using the Glauber-Mueller 
formula @ with Q sp . RdA shows a Cronin peak for k± in the range of Q s . We have used a prescription 
such that the region z ~ 1/A does not affect the integral. For that purpose, we make the replacement 
In j37\7 - * m (i3i\7 fl2 ) choosing a = 3 @]. A good approximation of the integral over b is to choose 
for Q s the average value (Q 2 (b)) = Q 2 min _ bias = (2/3)Q 2 c , in which case 

Actually, IjlUfl turns into an equality in the high k± region and is quite good in the region of the Cronin 
peak. In Fig. IHb) we see that the error is maximum when k < Q Sl reaching 10%. For a cylindrical 
nucleus there is no b dependence in the cross section and (|10J) turns into an equality. This tells us that 
the physics is the same whatever the geometry of the nucleus [H] • 

2.1.2 Effects of fragmentation on the Cronin peak 

To get the hadron cross-section we still have to convolute J7J with the proper fragmentation functions 

da dA^hx fdz,. „ 2 ,da dA ^ x . 
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Figure 1: (a) : The nuclear modification factor RdA without fragmentation functions (gluon production) 
for Q 2 C = 2 GeV 2 . (b) : Comparison between the average full min-bias calculation and the calculation 
with average min-bias Q s for Q 2 C = 2, 5 and 9 GeV 2 (Full, dashed and dotted). 



where Qf is a large scale of order k± taken between 2 and 8 GeV. We use the fragmentation 
functions of [2DI - I n Figl2 a ) we present the result for RdA- 

The fragmentation functions induce a flattening of the Cronin peak toward 1, with a shift toward the 
small fcj^'s. This feature can be understood by comparing R 9 dA without F.F's (fragmentation functions) 
and R'j A with F.F's, order by order in powers of Q 2 /k 2 (or the number of participants). Above the 
saturation scale, the following expansion has a meaning 

Ka = ^ + ^1 3 + - (12) 

In the Leading Log approximation 
modification factor by 



where we have defined for n ^ 2 

(z n )= [ dzD{z,Q 2 f )z n / [ dzD(z,Q 2 ). (14) 

Since z < 1 we always have (z n )/(z 2 ) < 1. The rescattering terms are then less important and R dA 
gets closer to 1. This is shown in Fig. [21 (a) : comparing with Fig. ^ (a) we see the dramatic effect of 
the fragmentation functions for Q 2 min _ bias = 1.3 GeV 2 . To get agreement with RHIC data we have to 
increase the value of the saturation scale Q s and take Ql. m i n -bias ^ 6 GeV 2 . This feature has already 
been mentioned in 0. In fact, Q 2 as defined in (JBJ is at most of the order of 2 GeV 2 for b = 0. The 



with respect to ln(fcj_/A), using ljll|l . we approximate the nuclear 
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Figure 2: (a) : Minimum bias RdA for charged hadron production for Q 2 S 
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GeV" (lowest to highest curve) . (b) : Minimum bias RdA for neutral pions production for Q 2 
GeV 2 . The points are representative data taken from [5] . 
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authors of [7] propose to enhance it by adding an additional momentum due to " non-perturbative" 
nuclear effects : 

QS-QS + kW, (15) 

with k 2 = 1 GeV 2 . This amounts to boosting Q 2 (b = 0) to ~ 9 GeV 2 . In fact, as shown in Fig. [21(a), 
we need even a larger saturation scale. On the other hand, the way RdA is normalized (by estimating 
the proton-proton cross-section using (JjJJ) is not fully convincing. The ratio Rcp does not suffer from 
the same uncertainty. It is shown in Fig. El for Q 2 C = 9 GeV 2 and Q 2 sC — 2 GeV 2 for central and semi- 
central collisions. Taking into account the experimental error-bars, the large value for Q s c is definitely 
preferred and in agreement with data. 

We should at this point remark that saturation physics in the present stage is in the situation where 
the leading order perturbative QCD description was, concerning large k± spectra in hadron-hadron col- 
lisions. The agreement with data could only be obtained by implementing an intrinsic non-perturbative 
transverse momentum for partons inside the hadron. The present state seems to be that, at next-to- 
leading order, the perturbative theory becomes predictive |^ |221 : a detailed NLO comparison with 
RHIC data for pp — > ir°X is presented in |23|. In this respect and in the leading order perturbation 
QCD framework, the authors of ref. [H] have used a traditional Glauber-Eikonal approach of sequen- 
tial multiple partonic scatterings with the implementation of a large intrinsic k± in parton distribution 
functions and they have obtained a good agreement with data for no production in dA collisions. 

Staying in the saturation physics framework, we may nonetheless try to go beyond the ad-hoc 
ansatz Ijl5|l and modify the saturation picture, which is exclusively based on hard multiple scatterings, 
by adding non-perturbative scatterings, when a parton (gluon) is passing through a nucleus. Following 
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Figure 3: Rcp for Q 2 sC = 9 GeV 2 (thick lines) and Q 2 sC — 2 GeV 2 (thin lines). Full lines correspond to 
central over peripheral collisions (full experimental dots). Dashed lines correspond to semi-central over 
peripheral collisions (empty experimental dots). Data from [2]. 



the Moliere scattering theory [21] extended to QCD [25], we define the probability distribution for the 
scattered parton by 

„„ , 1 da 1 , c —4- (1 - c)A 2 . , N 

v y crd 2 k 7r L (k 2 ) (k 2 +A 2 ) 2J v y 

with / d 2 kt/(k) = 1. 

A soft, non-perturbative, gaussian contribution is added to the hard screened (by mass A) gluon 
exchange term. Solving the kinetic master equation for the survival probability of the propagating 
parton, an effective scale (up to logarithm) is derived, 

Ql^(c^ + (l-c))Ql (17) 

Different from H15[) the non-perturbative part is added but weighted by a factor c < 1. 

In order to obtain an effective Q 2 S O(10 GeV 2 ) one has to add a significant non-perturbative part, 
e.g. for (k 2 ) ~ 0.5 GeV 2 , A = 0.2 GeV and c ~ 0.3 indeed 

Ql sa 2 GeV 2 -» Q 2 S w 9 GeV 2 . (18) 
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The underlying picture of dominating soft parton interactions with nucleons has strong implications, 
especially for the width of jet broadening, i.e. the resulting transport coefficient for cold matter q ~ 
Q 2 /(2R) becomes also rather large, namely q ~ 0.8 GeV 2 /fm (for gold)! But such a strong jet broadening 
has not been observed so far |27| . 



2.2 Forward rapidity and quantum evolution 
2.2.1 Hadron production cross-section 

As shown in |HJ^j, it is possible to rewrite eq. J3J under a fc^-factorized form which is then 
generalized to include the rapidity dependence. At leading twist, on defines the unintegrated gluon 
distribution for the nucleus and the deuteron, respectively, as 

Va(Y a , k, b) = — ^f-g / d 2 zWlN G (Y A , z, b)e- ik z , (19) 

and 

cp p (Y d , k, b - B) = J d 2 zV 2 z n G (Y d , z, b - B) e -* k - z , (20) 

with Ya = \b.{1/xa) = Y + t] and Yd = lxi(l/xd) = Y — r\ the rapidities of the gluons merging respectively 
from the nucleus and the deuteron and carrying the light cone momentum fractions xa and xa', rj is the 
rapidity of the produced gluon measured in the forward deuteron direction and Y = \n(^/~s/k±). One 
obtains the expression for the cross-section as 



da dA-, g x ^ 2(Xs x 
drjcPkcPb ~ C/7k 2 



i J d 2 B J d 2 qtp d (<l ~ k, Y - r),b - B)^(q, Y + v, b). (21) 



We will focus our analysis on the large forward rapidity region for the gluon where we can neglect the 
emission of additional gluons in the wave function of the proton. Actually, in this region the biggest 
value of Yd is reached at 77 = yielding Yd — Y which is not indeed large enough at the energies of 
RHIC for the evolution to take place. However, the opposite happens in the nucleus where Ya increases 
with 77 and we thus expect to probe the small x regime of the nucleus wave function K We have, taking 
(fid — 2(p p at leading twist without quantum evolution |18j . 

|d 2 bMk,y-^0,b) = ^A. (22) 

This allows us, assuming ipA smooth enough, to approximate l|21|) at leading log accuracy, when |k| ^> A, 

by 

da dA^ g x 2a s tp A (k,Y + r h b) f - // |k| . n „, - 

dfj^h=C^ k 2 J dB J A M*Y-^ ,B)d\ (23) 

Using the relation, 



we end up with 



xG(x d ,k 2 )= J d 2 B <p d (q,Y- V ~0,B)^, (24) 
da dA^ g x as (2 7r ) ¥ , A (k,y + r? ! b) 



dr]d 2 kd 2 b C F k 2 



-xG(x d ,k 2 ). (25) 



calculation of RdA has been recently performed 1201 in the framework of linearly factorized pQCD at NLO, taking 
into account nuclear leading twist shadowing for the partonic distributions in the nucleus. The observed suppression 
cannot be explained in this context. 
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Notice the collinear factorization of the gluon distribution in the proton. For hadron production one 
has to add the valence quark contribution which may be written [S] 

da dA-* q x a s (2Tr) <p A {k,Y + r},b) _ ______ l2 



xq v {x d X). (26) 



dr]d 2 kd 2 b N c k 2 
Convoluting with the fragmentation functions we get 

da dA ^ hx a s (2ir) ^ f 1 y A (k/z, Y + r, + In z, b) r 2 2 2 

where f Uj d{x,k 2 ) = (Cp/N c )xq U;l i(x,k 2 ) and / g (x,k 2 ) = xG(x,k 2 ) are the parton distributions inside 
the proton; D^/^Zjk) are F.F's of the parton i into hadron h, and zq — (k±/ y r s)e r > . Notice that this 
formula is identical to eq. (22) of I15j in the leading twist approximation. At forward rapidity, the 
longitudinal momentum fraction Xd carried by the parton inside the deuteron is of order of one. This 
implies that the dynamics of the deuteron is completely dominated by valence quarks. For numerical 
calculations we use the GRV LO parton distribution functions inside the proton from |28| , assuming that 
there is no significant difference between the proton and the neutron for charged hadron production. 

2.2.2 The unintegrated gluon distribution 

In the physics of saturation the relevant observable is the forward scattering amplitude of a quark- 
antiquark dipole off a target (a nucleus in our case). It enters several processes at high energy like DIS, 
photoproduction and hadron-hadron scattering. Its quantum evolution has recently been the object 
of many studies [El HI HE EH E2 ■ The Balitsky-Kovchegov (BK) equation [TH] valid in the 
large N c limit and in the mean field approximation provides a tool to study the rapidity behavior of 
the gluon distribution in a large momentum range including the saturation region where the effects 
of gluon recombination, taken into account in the non-linear term, become important, instead of the 
BFKL equation |34| which contains only gluon splitting (linear term) yielding a power growth of the 
gluon distribution with respect to the center of mass energy violating unitarity and the Froissart bound. 
In this work, as a simple test of the theory, we focus on the region above the saturation scale. The BK 
equation reads in momentum space 

d y N{k,y) = ^[ x (-^^)N(k,y)-N 2 (k,y)}, (28) 

where x{l) — 2i/j(1) — V'(t) — "0(1 ~ 7) is the BFKL kernel, and 

A>(k, b,y)= [ ^N(z, b, y)e lzk . (29) 



z 



At large y, a solution has been found for the BK equation above and not too far from the saturation scale 
in terms of travelling waves. It has a geometric scaling behavior in the variable L = ln(k 2 /Q 2 (b, y)) 
when y goes to infinity |3(JI I32| 

N(L, Y) = CL exp[- 7s i - f3(y)L 2 } . (30) 

C is an undetermined constant irrelevant in the present analysis of RdA and Rcp, 7s — 0.628 is the 
anomalous dimension of the BFKL dynamics in the geometric scaling region 123 Q 2 S S k± ^ 
Ql{Y) exp (1//3) and (3 = /3{y) = (2a X " (is)yy 1 ■ We used the recent fit to the HERA data performed 
in ESI where (3 = {2\ny)-\ A = 0.25, k = 9.9 and Q 2 s (y,0) = 3/2Q 2 smm ^ as (y) = 3/2A 1 /3( Xo / x )A 
GeV 2 with xq = 0, 67. 10 -4 and x = e~ v given here by x — k±/ \/se~ n . There is no straightforward way 
to link this asymptotic form to the one at mid-rapidity: this implies that there is an overall constant 



8 



coming from Q s which is negligible at very large rapidity and far from the saturation scale. We expect 
this constant to play a significant role at RHIC energies. Then, one can fix Q s as shown above, and put 
all the freedom into an additional constant (L — ► L + Lq). Thus, 

N(L, y) = C'{L + L ) exp[-( 7s + 2[3L Q )L - PL 2 }. (31) 

Making the same approximations as in the BK equation, one can write an expression for N G 6 

N G (z,y) = ^(2N(z,y)-N 2 (z,y)). (32) 

In the leading twist approximation where one can neglect non-linear terms, we find 

N G (z,y)~2N(z,y). (33) 

Both distributions N and ifiA are linked to the dipole scattering amplitude in (|29|l and (|19H . One can 
thus eliminate the latter yielding 

At very large y one can neglect the term which breaks scaling, namely (3 <C 1. So that the last expression 
reduces to the simple form exhibiting an exact scaling behavior 

2 

<Pa(L, y) cc (L + L )exp[-7 s L]. (35) 

7s 

We have fixed Lq such that tp A has a maximum when k 2 = Q 2 (h, y) corresponding to Lq = 3/7 s . 

This is the only free parameter of our calculation, it exhibits the uncertainty in the value of Q s . It turns 
out that Rcp is very sensitive to variations of Lq at energies of RHIC. For the numerical study we 
choose three different expressions for tpA, selecting various specific terms in <|31[) : 

i) The BFKL saturation-inspired form |33| . which violates scaling, derived from 

N(L, y) BF KL oc exp[- 7s L - f3L 2 }. (36) 

ii) The BK exact-scaling form (|35fl . 

iii) Finally the full expression derived from l|31|l . 



In Fig. 0] we show the comparison between the BFKL + saturation form H3(jfl and the BK 
parametrization 131|) . We expect this comparison to be valid at large enough k±. The agreement 
with data is quite good for the latter. With decreasing rapidity we would expect our formula to break 
down, nevertheless the global features of the data are reproduced even at 77 = 1. The exact scaling 
form, shown in Fig. |SJ is a slowly varying function of 77 and is too low to describe the data. However, 
if the picture is right, for increased rapidity (77 ~ 5 or 6), data points should match that shape. The 
fact that the saturation model has a semi-quantitative agreement with data for R G p is essentially due 
to the anomalous dimension since an approximate form is (when k± > Q s ) 

R CP ~ ^a#)V/ ~ {-&±Y»"-\ (37) 

coll part part 

At forward rapidity Jeff — 7s + fiiv) ^(^l/Qs) i s a decreasing function of 77 and an increasing function 
of k±. This allows us to understand the qualitative behavior shown by data and in particular the 
inversion of the centrality dependence compared to mid-rapidity where Rcp ^ 1 (Cronin enhancement) 
corresponding to "7e//" ^ 1 • At very large 77 the anomalous dimension stabilizes at Jeff = 7s, which 
could be tested at the LHC. 
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Figure 4: Rcp for the BK parametrization (thick lines) and the BFKL + saturation form (thin lines) 
at different rapidities r/ = 1,2.2 and 3.2. Full lines correspond to central over peripheral collisions 
(full experimental dots). Dashed lines correspond to semi-central over peripheral collisions (empty 
experimental dots ). Data from [2]. 
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Figure 5: Rcp for the exact scaling BK form at 77 = 
semi-central/peripheral. 

3 Conclusion and outlook 



3.2. Full line : central/peripheral, dashed line: 



To conclude we may reiterate that in our opinion, the saturation picture as probed by mid-rapidity data 
for the nuclear modification factor at RHIC energies cannot be seen as predictive. Indeed, the CGC at 
mid-rapidity is based on leading order pQCD calculations including multiple eikonal partonic rescatter- 
ings inside the nucleus, giving rise to a saturation scale Q 2 S at most of order of 2 GeV 2 and it turns out 
that this is not sufficient to describe RHIC data at mid-rapidity. It may be very interesting to relate 
the increase of the saturation scale needed to explain the dA data to the comparison of experimental 
measurements with theory for observables such as jet broadening in nucleus- nucleus (AA) collisions and 
also jet quenching, since the saturation scale determines their order of magnitude. 

At forward rapidity, a saturation inspired framework for quantum evolution predicts the suppres- 
sion as observed in data. This is a leading twist effect driven by the anomalous dimension which is 
the main source of the observed behavior. There is good evidence that the answer to the question: 
"has saturation been observed at RHIC ?" is positive and a definitive answer could be provided if one 
measures, in a more precise way, the anomalous dimension as predicted by the theory. To do so, one 
needs to go to larger rapidities (energies) at LHC. 

As one more step to improve further the status of the present theoretical description let us remark 
that initial state suppression effects may be also present in AA collisions at large rj and large transverse 
momentum k±. An interesting quantity to be measured is the double ratio R(r]i, T]2,k±) with 771 > 772 > 
0, defined by 

R(Vl,V2,k_L = - — j— , (38) 

Rcp (?72,k_l) 
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as a function of k±. For 771 = 2.2 and 772 = 0, this ratio is measured for AA by the BRAHMS Collabo- 
ration [371; there is evidence of suppression becoming more important at forward rapidities, as observed 
in dA measurements. 

Because of the final state suppression in AA collisions one may argue that this ratio is bounded by 

R(r)x,r)2,k±) < R(r)ur]2,kj_)i n iuai < 1, (39) 

where i?(?7i, 772, k±)i n i t i a i for AA collisions may be estimated for Rep-initial in a way done for dA colli- 
sions (e.g. by assuming k± — factorization). We expect it to be quantitatively similar to the dA case (e.g. 
Figs. 0Jb) and (c)), which is well understood in the saturation picture for 771 > 772 > 0. The inequality 
(l3~9"|) is based on the assumption that Rcp in AA collisions may be factorized as Rcp-initiaiQ(k±,ri): 
the quenching factor takes into account the (radiative) final state interactions It depends, at fixed 
fcj_ on the pathlength of the parton propagating in a dense medium, such that Q is decreasing with 
increasing pathlength, i.e. it amounts to more quenching. Assuming, on geometrical grounds, a longer 
path for 771 than for 772, 771 > 772, implies Q(fc±, 77i)/Q(fc^, 772) < 1, and l|39(l follows. 
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